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We argue that a large negative running spectral index, if confirmed, might suggest that there 
are abundant structures in the inflaton potential, which result in a fairly large (both positive and 
negative) running of the spectral index at all scales. It is shown that the center value of the 
running spectral index suggested by the recent CMB data can be easily explained by an inflaton 
potential with superimposed periodic oscillations. In contrast to cases with constant running, the 
perturbation spectrum is enhanced at small scales, due to the repeated modulations. We mention 
that such features at small scales may be seen by 21 cm observations in the future. 



I. INTRODUCTION 

Cosmic inflation [IH3] sets the initial conditions for the 
subsequent Hot Big Bang cosmology, and also generates 
the primordial curvature perturbation of the universe. 
The nature of this primordial density perturbation is 
studied through various cosmological observations, and 
especially CMB data can be well fitted by a power-law 
primordial spectrum, which is a typical prediction of the 
standard slow-roll inflation. 

Interestingly, however, when relaxing the assumption 
of a power-law spectrum, latest CMB observations show a 
slight preference for a fairly large negative running spec- 
tral index. Ignoring tensor modes, the WMAP 7-yr data 
combined with ACT 2008 data give the bounds rig = 
1.032±0.039 anddns/dln/c = -0.034±0.018 (68% CL) at 
the pivot scale = 0.002 Mpc~^ (4j. Even with the ten- 
sor modes, similar bounds are obtained using the WMAP 
data only [5] . In the usual slow-roll inflation model with a 
featureless potential, the running spectral index is second 
order in the slow-roll parameters, and so it is expected to 
be of O(10~^). Thus, although the evidence for the run- 
ning spectral index is not conclusive, it is worth studying 
its implications for the inflation models. 

The center value of the running spectral index, if taken 
at face value, suggests either the inflaton potential with a 
step-like feature (see e.g. [6HTT] and references therein), 
or multiple inflation jl2|, I13j . since otherwise inflation 
would terminate within about 30 e- foldings [Hj . In both 
cases, fine-tuning is needed to ensure that the CMB scale 
exits the horizon at a specific time when the inflaton 
comes close to the feature in its potential, or at about 
10 e-foldings before the end of the first inflation. Such 
fine-tuning seems a puzzle and calls for explanation. 

We argue in this paper that, if there is indeed large neg- 
ative running spectral index, this might suggest that the 
inflaton potential has structures everywhere in the po- 
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tential and such a potential may be common in the land- 
scape of inflation models. Such repeated modulations 
to the inflaton potential can account for a perturbation 
spectrum which is blue tilted with negative running at 
the CMB scale, if one period of the modulations is large 
enough to encompass the observed scales. Inflaton po- 
tentials with modulations can lead to repeated features 
in the perturbation spectrum, in contrast to the case of 
single-stepped potentials. See also [T3] for related work 
in the context of Natural Inflation models (THl [13 with 
extra dimensions [HI [19] . 

Modulations to the inflaton potential can arise due 
to the microphysics governing inflation. Especially for 
large-fleld models, one expects that the super-Planckian 
held ranges may contain many features whose scale is 
typically related to that of the underlying microphysics. 
Explicit constructions of large-field models in string the- 
ory [2nH2^ actually show that periodic modulations of 
the potential can show up, whose period is determined 
by e.g., the size of the internal compactified space. (See 
also ^23', '2T for discussions on small-scale features along 
the infiaton field range in the cosmic landscape.) 

To be explicit, we consider large-field infiation from 
power-law potentials with superimposed periodic modu- 
lations, and show that the center value of the running 
spectral index can be realized. What differs from the 
previous studies using single-stepped potentials is that 
the running spectral index is not localized, but it is neg- 
ative at the entire CMB scales and becomes positive at 
small scales. Such characteristic behavior of the density 
perturbation may be seen by the future 21 cm observa- 
tions [IB] . The potential detection of such oscillatory 
feature in the power spectrum by future observations was 
studied in Ref. [27j, with the main focus on oscillations 
with (much) shorter period. 

The rest of the paper is organized as follows. In Sec. |ll] 
we briefly explain our idea of generating the running 
spectral index. Then we estimate the power spectrum of 
the density perturbation based on inflation models with 
a linear, and then a generalized polynomial potential su- 
perimposed with periodic modulations, both analytically 
and numerically in Sec. |III| and Sec. |IV| The last section 
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is devoted to discussions and conclusions. 



II. BASIC IDEA 

Let us first give our basic idea in this section. In this 
paper we use the units c = h — Mp — 1. 

The fact that the CMB data can be well fitted by a 
Gaussian power-law primordial spectrum indicates that 
the standard slow-roll inflation paradigm is a valid de- 
scription of the evolution of the early universe. However, 
if one looks at the details of the inflaton potential, there 
might be non-trivial structures which do not change the 
overall behavior of the inflaton, but nonetheless affect the 
tilt and/or the running spectral index in a non- negligible 
way. 

If the running of the spectral index is large at the CMB 
scales, it may suggest that the inflaton potential has step- 
like structures everywhere, since severe flne-tuning would 
be needed otherwise. So we consider an inflaton potential 
V{(jji) with modulations such as shown in Fig. [l] We 
decompose the potential as 




V{^)^Voi(P)+Vrr,odW, 



(1) 



FIG. 1: A sketch of the inflaton potential with superimposed 
periodic oscillations. The dashed line is a smooth potential 
without modulations. The modulations are amplified for visu- 
alization purposes. The resulting curvature perturbation and 
its tilt at scales exiting the horizon as the inflaton approaches 
the markers are shown in Fig. [2] and Fig. |3] The observed 
CMB scales are considered to lie within a half oscillation pe- 
riod. 



where the flrst term is a smooth potential sourcing a 
(nearly) constant spectral tilt, and the second term rep- 
resents the modulations. We assume that the overall in- 
flaton dynamics (such as the Hubble parameter and the 
inflaton velocity) is not altered drastically by the modu- 
lations, i.e.. 
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We further suppose that the slow-roll approximations 
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are valid. Then, in order for the running of the spectral 
index to be large enough so that the spectral tilt switches 
between blue and red within the observed scales, we re- 
quire 
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Thus, we are invoking a hierarchy 
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We note that the conditions ([5]), (|6]), and ([7| need not 
hold all along the inflaton potential, but are only required 
to be satisfled occasionally during inflation for large run- 
ning of the spectral index to be produced. 

An important assumption here is that the effect of the 
modulations Vmod (and its derivatives) on the inflaton 
dynamics is sufficiently small when averaged over a suf- 
ficiently long time or large field space. One example for 



such modulations is a sine function, and we will see in the 
following sections that the above conditions can be easily 
satisfied in this case. We emphasize here that the effect 
of modulations can be missed in such an analysis that the 
inflaton potential is expanded in terms of a power of up 
to a flnite order |14 . On the other hand, as we shall see, 
the spectral index and its running are easily modifled by 
the modulations, while not affecting the overall inflaton 
dynamics. 

The spectrum of the curvature perturbation can be 
computed, 
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which should be estimated at the time of horizon exit k = 
aH. Note that, since (j> is determined by V , the 
main contribution to the power spectrum arises from the 
smooth part Vq, and that the V^^^^ results in small os- 
cillations of the perturbation spectrum as illustrated in 
Fig. [2j The oscillatory feature of the perturbation spec- 
trum arises mainly from </) (i.e. V^^^), and the modu- 
lations of H (i.e. Vmod) is subdominant in the following 
analysis. The spectrum is blue (red) for scales exiting 
the horizon when the inflaton is rolling along the poten- 
tial whose curvature is positive (negative). To be precise, 
since the dropping of the Hubble parameter reddens the 
perturbation spectrum, the resulting spectrum becomes 
blue tilted when the deceleration of the inflaton (i.e. the 
positive curvature of the potential) is large enough. 

Before computing the tilt of the spectrum, let us pause 
and comment on the validity of the slow-roll approxima- 
tions. Expressions for the spectral index and its running 
(and its running, and so on) derived using the slow-roll 
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FIG. 2: The curvature perturbation spectrum generated from 
the inflaton potential in Fig. [l] 
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FIG. 3: Spectral index and its running from the inflaton po- 
tential in Fig. [l] 



approximations contain errors of order of the approxi- 
mate results multiplied by 
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where 
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In the final expressions, we have kept contributions from 
the smooth potential Vq which are denoted by parame- 
ters with the subscript 0, and also the leading oscillatory 
terms due to the modulations Vmod, by using the hier- 
archy Q and assuming that the inequality ([3| is well 
satisfied. When the oscillatory terms are dominant, the 
tilt of the perturbation spectrum repeatedly switches be- 
tween red and blue, hence the spectral index — 1 as 
well as the running dris/dlnk change between positive 
and negative values. This results in a clockwise motion 
in the (us — l) - dus/dhik plane as inflation proceeds, 
as illustrated in Fig. [3j Given that the pivot scale ex- 
ited the horizon when the inflaton was around the filled 
circle in Fig. [T] in one of the repeated modulations, a 
blue tilted perturbation spectrum with negative running 
is obtained. 

The hierarchy ([t]) implies that even if the modulations 
are so small that their effects on the Hubble parameter 
and the inflaton velocity are subleading, they can still sig- 
nificantly contribute to the spectral tilt of the perturba- 
tion spectrum. However, one should also note that such 
hierarchy can lead to the breakdown of the conditions ^ 
and ^ at other places along the inflaton potential, thus 
considerably affecting the inflaton dynamics in the end. 
We will see in the following sections using explicit exam- 
ples that oscillatory modulations can satisfy the above 
conditions all through (or at least for a large fraction of) 
the inflationary period. 
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III. LINEAR POTENTIAL WITH PERIODIC 
OSCILLATIONS 



Therefore, as long as this is smaller than unity, results ob- 
tained using the slow-roll approximations should be good 
estimates. Such condition is satisfled for the explicit cases 
we study later on, since we require the observed scales to 
be generated within a half oscillation period. For de- 
tailed analyses of inflation and primordial perturbations 
beyond the slow-roll approximations, see e.g. [28H35] . 

Then the spectral index and its running can be esti- 
mated as follows, 
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Now let us consider an explicit inflation model. As 
a simple example, we consider curvature perturbations 
from a linear inflaton potential |21l 136) with superim- 
posed periodic oscillations [22j|37] (cf. Fig. [T]): 



V{(j)) = + cos 



f 



(13) 



where A, A, /, and 9 are constants. The period of the 
oscillations is determined by /, which we assume to be 
sub-Planckian, i.e. / < 1. This is often the case in ex- 
plicit microscopic constructions, and is also compatible 
with our requirement that half of an oscillatory period 
roughly corresponds to the time duration when the ob- 
served scales exited the horizon. We also assume V > 0, 
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namely, 



A/ 



< 1, 



(14) 



which guarantees that the inflaton is not trapped by the 
modulations during inflation. This condition is equiva- 
lent to (|3|. 

Assuming that the slow-roll approximations Q are 
valid, the spectrum of the curvature perturbation can 
be computed, 



A4 



(15) 



which should be estimated at the time of horizon exit. 
Note that in obtaining the right hand side we have omit- 
ted modulations of H (i.e. V) which is subdominant 
due to ( 14 1 and f <^ (f> (since inflation occurs for super- 



Planckian tf>). 

Let us now check whether we can use expressions for 
the spectral tilt derived using the slow-roll approxima- 
tions. (See also discussions below Upon its compu- 
tations, one finds 



V d 
Vdch 



< 



1 



(16) 



where the / in the denominator of the right hand side 
comes from differentiating the oscillatory part. As we will 
soon see, in our case {4>f)~^ is smaller than unity, at least 
until when the CMB scale exits the horizon, due to the 
requirement that the period of the oscillations be large 
enough to encompass the observed scales. Therefore we 
trust results derived using the slow-roll approximations. 

Then the spectral index and its running can be esti- 
mated as follows, 
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In the final expressions (17) and (18), we have only kept 



the leading oscillatory terms. One clearly sees that the 
spectral index and its running oscillate with growing os- 
cillation amplitudes as inflation proceeds and 4> becomes 
smaller. (The relative amplitude compared to the non- 
oscillatory part becomes smaller.) Note that the power 
of /(< 1) appearing in Eqs. ( 17 ) and ( 18 ) are larger than 
that of Eq. (15), and that is why the spectral index and 



its running can be significantly modified without affect- 
ing the power spectrum too much. One can easily check 
that the conditions ([2]) - ([6| can be satisfied for a certain 
choice of /. 

Let us now set the parameters of the potential so that 
the spectral index and its running at the pivot scale re- 



alizes the central values of the bounds from WMAP 7- 
yr [5]: 



n, = 1.076 ±0.065, 
dris/dlnk = -0.048 ±0.029, 



(19) 



(both at 68% CL) when both the running and the tensor 
mode perturbations are allowed to vary, in addition to 
the spectral index. (The bound on the tensor-to-scalar 
ratio here is r < 0.49 (95% CL).) Note that the pre- 
dicted power spectrum is not exactly the one used in the 
analysis by the WMAP team. However, since the pe- 
riod of oscillations is sufficiently long compared to the 
observed CMB scales, such a difference will not signifi- 
cantly change the likelihood distribution. For simplicity 
we use the above values in the following analysis. 

We set the pivot scale to have exited the horizon when 
the inflaton field was at around (/jq — 10, which corre- 
sponds to about 50 e-folds before inflation ended. For 
a linear potential (without modulations), the number of 
e-folds obtained within the field range A0 starting from 
00 is AA/" ~ 4>o^4>- Here we set / by choosing the half 
period of oscillation from the pivot scale to correspond 
to about In^CMB ^ 7 5 e-foldings with ?Cmb _ 2OOO, i.e. 
tt/ = AJ\f /20O = 0.75. (Note that requiring the observed 
scales to lie within one oscillation period, i.e. AAf > 10, 
leads to {4>of)~^ < 1 which validates the analytic ex- 
pressions derived using slow-roll approximations, as was 
discussed below ^ and (16).^) Then the remaining pa- 
rameters A, A, are determined so that the perturbation 
spectrum at the pivot scale realizes the COBE normal- 
ization [38j = 2.43 X 10~^, as well as the central values 
of We choose A « 1.99 x IQ-io, A « 1.77 x 10-3, 
6 « —2.03, and the inflaton field value when the pivot 
scale exits the horizon as </> « 9.98. We note that the 
modulations' amplitude A is of order the GUT scale. 

With these parameter sets, we have numerically solved 
the equation of motion of the inflaton and the Friedmann 
equation, and computed the curvature perturbation spec- 
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which is shown in Fi; 
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cillatory features are clearly seen in the spectrum 
comparison, we also plotted a spectrum with constant 
running a, 



Pcik)^Pc{ko) 



ko 



(20) 



extrapolated from values at the pivot scale ko — 
0.002 Mpc^^. The CMB observations probe scales from 
k ~ lO""* Mpc""'^ down to ^ 0.1 Mpc"^, and within this 
range these two lines are almost degenerate. The differ- 
ence becomes evident only at small scales fc > 10 Mpc"^, 
which may be probed by the future proposal of the 21 cm 



^ The condition (</>/) ^ < 1 is satisfied until about 10 e-foldings 
before the end of inflation. 
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FIG. 4: Curvature perturbation spectrum (black solid line) 
from the inflaton potential (13 1, realizing central values of 
(19 1 at the pivot scale fco — 0.002 Mpc~^ which is denoted by 
the filled circle. Spectrum with constant running (201 extrap- 



olated from values at the pivot scale is shown as the black 
dotted line. Spectrum from a linear potential without any 
modulation is shown as the black dashed line for comparison. 



observations [551 US] ■ We also plot the varying spectral 
index and its running in Fig. [5] For each oscillation in 
the inflaton potential, the trajectory in the phase space 
passes through the WMAP constrained region. As in- 
flation proceeds, i.e. towards smaller scales, the radius 
of the orbit increases and also its center shifts towards 
smaller n,. 



IV. POWER-LAW POTENTIAL WITH 
PERIODIC OSCILLATIONS 

All the previous discussions can be straightforwardly 
extended to large-field inflation with general power-law 
potential (see e.g. [T3 EHl SI] for models which exhibit 
quadratic inflaton potentials with superimposed periodic 
oscillations). 




FIG. 5: Variation of the spectral index ris and its run- 
ning dus/dlnk as inflation proceeds along the modulated po- 
tential (131. The plot starts from 70 e-folds before the end 
of inflation (end point of the innermost orbit), up to ~ 10 e- 
folds before the end of inflation where most part of the orbit 
goes beyond the plotted region. Filled circle corresponds to 
the pivot scale (about 50 e-folds before the end of inflation). 
The scale which exited the horizon 15 e-folds before (after) 
the pivot scale is indicated by an open (filled) star. Super- 
imposed contours are constraints from 5-year WMAP (blue) 
and from 5-year WMAP-hBAO-f SN (red) (68% and 95% CL), 
when tensor moder perturbations, spectral tilt, and running 
index are allowed to vary |39) . 



ilar fashion. 
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Viq^) = A(/)" + A^ cos ( ^ + 



(21) 



The condition equivalent to (14 1, which guarantees V > 
is now 



A4 



< 1, 



(22) 



but one should note that for n > 1, this may break down 
as (j) becomes smaller and the inflaton can be trapped in 
a local minimum due to the modulations. 

The approximate expressions for the perturbation 
spectrum and its spectral tilt can be derived in a sim- 
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One can see that the spreading of the oscillation ampli- 
tudes of ris— I and dus/dln k as inflation proceeds is more 
drastic for larger n, since oscillations with constant am- 
plitude as in (21) become relatively more significant at 
small <j). For comparison, we have plotted oscillating tra- 
jectories in the (ng — 1) - dus/dhik plane for the linear 
{n — 1) and quadratic (n = 2) potentials with modula- 
tions in Fig. [6j The parameters of the potentials are cho- 
sen in order to realize the central values of ( 19 1 and the 
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FIG. 6: Variation of the spectral index Us and its run- 
ning dus/lnk for linear (solid) and quadratic (dashed) po- 
tentials with modulations. Filled circle indicates when the 
pivot scale exited the horizon. Both trajectories are plotted 
from 15 e-folds before the pivot scale exited the horizon (end 
points of the innermost orbits), up to 20 e-folds after the pivot 
scale (end points of the outermost orbits). 



COBE normalization at the pivot scale. The case with 
the linear potential was discussed in the previous section. 
For the quadratic potential, the pivot scale is set to have 
exited the horizon at around (f>o — •\/200, and / is chosen 
to realize a half oscillation period lasting for about 7.5 
e-foldings from the pivot scale, i.e. tt/ = AA/'/0o ~ 1.06. 
The remaining parameters are chosen as A ~ 1.95x 10^^^, 
A w 2.52 X 10-3, e sa -2.06, and the inflaton field value 
when the pivot scale exits the horizon as « 14.1.^ 

Normalized at the pivot scale, one sees that the tra- 
jectories for the quadratic potential, compared to that of 
the linear case, starts from an orbit with smaller radius 
and grows to a larger orbit towards the end of inflation. 



V. DISCUSSIONS AND CONCLUSIONS 

Compared to cases with constant running of the spec- 
tral index, the perturbation spectrum from periodic mod- 



^ For this parameter set, the breakdown of the condition l |22| hap- 
pens at around the end of large-field infiation, hence we disregard 
the trapping of the inflaton in a local minimum due to the mod- 
ulations. However, modulations do flatten the potential towards 
small and extend the period of infiation: The number of e- 
folds obtained from when ipo ft: \/200 until the end of inflation is 
about 50 for a quadratic potential without modulations, but in 
this case about 57. Such increase/decrease of e-foldings due to 
modulations may have interesting consequences, see discussions 
in Sec. |V] 



ulations which we have studied is enhanced at small 
scales, thus does not delay structure formation [42]. 
Throughout this work we have assumed that the slow-roll 
approximations are valid during inflation. If we consider 
an extreme case where the slow-roll approximations are 
violated temporally during inflation, a perturbation spec- 
trum with large oscillations at the corresponding scales 
is realized. Such a spectrum may be used to account for 
the anomalous suppression at the quadrupole (see [i^Hi] 
and references therein) as well as the missing satellite 
and/or core-cusp problems [151 

We focused on a case that a sizable running of the 
spectral index is generated due to the modulations. We 
may use such a structure to modify the prediction on 
the spectral index to make an inflation model consistent 
with observations. One example is chaotic inflation with 
a quartic coupling, which is strongly disfavored by the 
WMAP and other observations because the tensor-to- 
scalar ratio is too large for the predicted spectral index. 
The tension can be relaxed if modulations terminate in- 
flation earlier compared to the case without modulations. 
Then the inflaton fleld value when the CMB scale ex- 
ited the horizon is corrected, modifying predictions on 
the spectral index and the amplitude of the tensor mode 
perturbations. Another possibility is that the modula- 
tions directly increase the spectral index. However for 
this to happen, the amplitude of the modulations should 
decrease towards smaller values of the inflaton, otherwise 
the inflaton would get trapped by one of the modulations 
as inflation proceeds. 

In this paper, we have explored the possibility of pe- 
riodic modulations to the inflaton potential giving rise 
to large running of the spectral index at all scales. Due 
to the oscillatory modulations, the spectral index — 1 
as well as its running dus/dhik switch between positive 
and negative values. Such behavior results in distinct 
features in the perturbation spectrum at small scales, 
which may be probed by future 21 cm observations. As 
explicit examples, we have studied inflation models with 
power-law potentials with superimposed periodic oscil- 
lations. We showed that the modulations can produce 
a running spectral index, while having small impact on 
the overall inflaton dynamics. In order for the cosmolog- 
ical observables to lie at the central values of the WMAP 
bounds, modulations with amplitudes of order the GUT 
scale were considered. Especially in large-fleld models, as 
the inflaton travels super-Planckian field ranges, modula- 
tions in the potential can show up with scales determined 
by the microphysics governing inflation. Concrete con- 
structions of such cases in the context of string theory 
were performed in |211 122] , which realized linear inflaton 
potentials with superimposed oscillations. 

Though we focused on modulations with a large os- 
cillation period encompassing the entire CMB scale, 
cases with small periods can also have interesting conse- 
quences, such as oscillations of the perturbation spectrum 
within the observed scales [571 SOI HZ] • Very small oscilla- 
tions whose period is sufficiently shorter than one e-fold 
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during inflation, can further lead to resonance effects in 
the power spectrum as well as in the non-Gaussianity, as 
was studied in [21 [371 Si US] . 

Large running of the spectral index, if confirmed, may 
suggest the presence of abundant features in the inflaton 
potential, indicating connections with the microphysics. 
It would be interesting to examine what we can say about 
the underlying physics if future observations detect com- 
binations of signals, such as tensor modes and a running 
spectral index. 
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